Towards a novel wave-extraction method for numerical relativity. 
I. Foundations and initial-value formulation 
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The Teukolsky formahsm of black hole perturbation theory describes weak gravitational radiation 
generated by a mildly dynamical hole near equilibrium. A particular null tetrad of the background 
Kerr geometry, due to Kinnersley, plays a singularly important role within this formalism. In 
order to apply the rich physical intuition of Teukolsky's approach to the results of fully non-linear 
numerical simulations, one must approximate this Kinnersley tetrad using raw numerical data, 
with no a priori knowledge of a background. This paper addresses this issue by identifying the 
directions of the tetrad fields in a quasi-Kinnersley frame. This frame provides a unique, analytic 
extension of Kinnersley's definition for the Kerr geometry to a much broader class of space-times 
including not only arbitrary perturbations, but also many examples which differ non-perturbatively 
from Kerr. This paper establishes concrete limits delineating this class and outlines a scheme to 
calculate the quasi-Kinnersley frame in numerical codes based on the initial-value formulation of 
geometrody namics . 
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I. INTRODUCTION 

After an early history marked by vigorous debate (see, 
e.g., P]), the existence of gravitational radiation has be- 
come accepted as a hallmark prediction of Einstein's 
general theory of relativity. Despite indirect evidence 
supporting this prediction, however, gravitational waves 
have yet to be observed directly. Indeed, rough, order- 
of-magnitude calculations indicate a typical passing wave 
will produce only tiny material strains of order 10"^"'^. 
A new generation of laser-interferometric gravitational 
wave observatories have recently been commissioned, and 
several are already operational, which aim to observe 
such strains and thereby detect gravitational waves 
incident on Earth. These experiments are necessarily ex- 
tremely delicate and, as a result, both initial detection 
and the long-term goal of extracting new information 
about the distant sources of a particular gravitational 
wave must be aided by detailed theoretical predictions. 

A community of theorists have turned to numerical 
simulations of the full, non-linear Einstein equations to 
identify the characteristic features of the gravitational ra- 
diation generated by various sources. An impressive ar- 
ray of techniques have developed within this young field 
of numerical relativity which aim to provide the ongoing 
experimental effort with accurate predictions. However, 
the program faces a number of challenges, from founda- 
tional questions to issues of implementation and interpre- 
tation of results. Here, we consider one issue of interpre- 
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tation. We ask how, exactly, the well-known theoretical 
description of gravitational radiation outside quiescent 
black holes due to Teukolsky ^ might be applied to the 
results of numerical simulations. 

The notion of gravitational radiation in general rela- 
tivity does not have universal meaning. Rather, a proper 
definition can be given only in regions of space-time, ra- 
diation zones, whose geometry is characterized by two 
distinct length scales one describing an "average" 
radius of curvature and a second, much shorter scale cor- 
responding to the wavelength of gravitational waves. Be- 
cause these two scales can be distinguished only when 
the waves contribute small corrections to the "average" 
curvature, many analyses of gravitational radiation are 
founded on perturbation theory. Examples include not 
only the standard analysis in linearized gravity, but the 
Regge-Wheeler and Zerilli approaches, as well as 
Teukolsky's. Even the asymptotic formulation of the 
gravitational radiation , which applies to quite general 
space-times, operates by studying differences between the 
physical metric and a fixed asymptotic metric near con- 
formal infinity. In numerical relativity, however, all of 
these analyses are difficult to implement. The perturba- 
tion approaches are complicated because no background 
metric on space-time is known a priori, while the asymp- 
totic approach needs access to field values far outside the 
computational domain. 

This paper focusses on Teukolsky's perturbative ap- 
proach since it describes radiation fields near a rotat- 
ing black hole, the expected end-state of many physical 
processes producing relatively strong gravitational wave 
signals. We first identify the essential elements of the 
Teukolsky formalism which derive from its fixed back- 
ground geometry. In particular, we are interested in the 
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preferred tetrad of basis vectors on space-time which un- 
derhes the definition of Teukolsky's fields. Although this 
tetrad arises as a perturbation of a canonical, Kinnersley 
tetrad ^ for the background Kerr geometry, we show 
one can approximate it using only the physical metric, 
eliminating any reference to the background. The con- 
struction of this approximate Kinnersley tetrad occurs 
naturally in two stages. The first, which is completed 
in this paper, fixes the directions of the desired tetrad 
vectors. Because the final results of our analysis make 
no mention of a background geometry, they may be ap- 
plied unambiguously to a broad class of space-times. We 
give concrete criteria characterizing this class, and find 
it includes many examples which differ non-perturbatively 
from the Kerr geometry. In particular, when a numerical 
simulation produces a single black hole settling down to 
equilibrium, this first stage should be viable even at rel- 
atively early times after the hole forms. The problem of 
the second stage, fixing the scalings of the tetrad vectors 
along these preferred directions, is described in some de- 
tail here but not solved. We plan to present a solution 
in a future paper. However, even the first stage alone 
provides significant gauge-invariant information partially 
characterizing the gravitational field encoded in the nu- 
merical variables. An earlier Letter presented one 
scheme to extract such information, and this paper can 
be regarded as an elaboration and expansion of that pre- 
sentation. 

Like the Teukolsky formalism itself, several of this pa- 
per's results are most simply expressed in the Newman- 
Penrose approach |0 to space-time geometry. However, 
because many existing numerical codes do not currently 
implement this approach, we strive to present our fi- 
nal results in a form which does not refer explicitly to 
Newman-Penrose quantities. Rather, we seek a form 
which would allow the Teukolsky fields to be deduced 
ab initio from the physical metric. There is a price to 
be paid for this, as certain results appear more compli- 
cated in this language. A complementary, and slightly 
more comprehensive, analysis based on the Newman- 
Penrose approach is given in a companion paper [lll| . 
This complementary analysis operates by picking an ar- 
bitrary tetrad on a generic space-time and showing how 
to transform it to yield one of interest in the Teukolsky 
formalism. We present the two analyses separately to 
maximize the clarity of each. 

The following is an outline of this paper. First, we re- 
view elements of the Teukolsky formalism from the tradi- 
tional point of view, where a preferred background Kerr 
geometry is known a priori^ and formulate in detail the 
problem to be solved. Second, we introduce the notion 
of a transverse frame. Several such frames exist in a 
general space-time, all are calculable from the physical 
geometry, and one is the quasi-Kinnersley frame of inter- 
est in the Teukolsky formalism. Finally, to help deploy 
these results for numerical relativity, we express these re- 
sults using only initial-value data on an arbitrary spatial 
slice. We conclude with several comments. Finally, an 



Appendix elaborates one of the central calculations of the 
paper. 



II. THE TEUKOLSKY FORMALISM AND 
NUMERICAL RELATIVITY 

The observable signature of a non-trivial gravitational 
field in general relativity is space-time curvature. How- 
ever, because no coordinate system is preferred over any 
other, the coordinate components of the curvature ten- 
sor do not generally have any invariant physical mean- 
ing. One must be careful to distinguish, for example, 
between actual gravitational radiation and oscillations in 
those components caused by a peculiar choice of coor- 
dinates. The Teukolsky formalism of first-order pertur- 
bation theory Q solves this problem using a family of 
gauge-invariant^ scalar variables originally proposed by 
Newman and Penrose (lOj . 

The Newman-Penrose approach introduces an initially 
arbitrary null tetrad n", m", to'*) on space-time. The 
four vector fields here are all null, with the first pair 
real and the second both complex and conjugate to one 
another. Their only non-vanishing inner products are 

rna = -l and rha = 1. (1) 

It is convenient to restrict attention to oriented null 
tetrads satisfying 

24z £[a m mc fhd] = £abcd, (2) 

where eabcd denotes the usual space-time volume element. 
This restriction can always be made without loss of gen- 
erality in cases of physical interest. The independent 
components of the Weyl curvature tensor in such a non- 
coordinate basis are encoded in five complex Weyl scalars 

^2-=Cabcd'i''ni^m^n'' „a b ^d (3) 

,T, ^ -a b-c d ^3--(^abcdi U m n 

These scalar quantities are naturally independent of any 
space-time coordinate system, but do clearly depend 
on the choice of null tetrad. Accordingly, the '^n can 
acquire direct (observer-independent) physical meaning 
only when a physically preferred null tetrad can be found. 

The Teukolsky formalism picks a particular null tetrad 
for the Newman-Penrose scalars by exploiting special 
features of its given Kerr background. Specifically, any 
Kerr geometry admits exactly two repeated principal 
null directions^ {^°} and {n"}. In any null tetrad 



^ Throughout this section, "gauge" refers to the first-order diffco- 

morphisms of linearized gravity. 
^ The braces {•} emphasize that only the directions of these vectors 

are determined. 
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(^°, n", m°, m") whose real null vectors lie in these di- 
rections, all the Newman-Penrose curvature components 
except \l/2 vanish, guaranteeing, as we discuss further be- 
low, the gauge-invariance of first-order perturbations in 
and 5*4. It is therefore quite natural to choose the 
directions of the tetrad fields in this way, but one must 
proceed to fix their scalings along these directions as well. 
A technique to do this in an exact Kerr, or generally in 
any Petrov type D, geometry was first proposed by Kin- 
nersley . To explain this technique, it is useful to recall 
the generic ambiguity in the null tetrad. 

At a given point of a general space-time, one can trans- 
form any given oriented null tetrad {£"■ , n'^ , m"- , fh°-) to 
any other using elementary transformations of three ba- 
sic types. These include two families of elementary null 
rotations L{a) and N{b), with a and b arbitrary complex 
parameters, which preserve either i°- or n", respectively, 
while changing the directions of all three of the other 
tetrad elements. A third family of spin-boost transfor- 
mations S(c), with c an arbitrary complex parameter, 
mutually scales the real null vectors and n°, leaving 
their directions unchanged, and rotates the complex ele- 
ments m" and m" by complementary phases. In addition 
to these, we highlight a single fourth transformation, the 
exchange operation E, which interchanges the real and 
complex tetrad elements in pairs. The detailed formulae 
describing these transformations are unimportant here. 
They are standard and can be found, for instance, in 
[Til Il3l|. We do, however, note the effects of these trans- 
formations on the Weyl scalars 
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where (^) denote the usual binomial coefficients. 

In the Kerr geometry, fixing £° and to lie in the 
two repeated principal null directions eliminates the null 
rotations completely, but does not restrict the spin-boost 
and exchange freedom in any way. To proceed, Kinners- 
ley implicitly breaks the exchange symmetry by choosing 

to point outward, toward infinity, in the exterior of the 
hole and then sets 

e:=ir(^''V,nb-m^V,mb) = 0. (5) 

The quantity e is one of the standard Ricci rotation co- 
efficients used in the Newman-Penrose approach. Set- 
ting it to zero does not eliminate the residual spin-boost 
ambiguity entirely, however, as one remains free to do 
spin-boost transformations with parameters c satisfying 

rv,c = o. (6) 

Note that, since the right side vanishes, this equation 
is independent of the scaling of the vector field and 
implies only that c must be constant along each of its 



integral curves. The set of such curves rule the exterior 
region of the black hole and each intersects a unique point 
of future null infinity. As a result, adapting the scaling 
of the tetrad near infinity to stationary observers there 
uniquely determines the Kinnersley tetrad everywhere on 
Kerr space-time. That is, fixing the scalings of and 
such that their spatial projections relative to a stationary 
observer have equal magnitude near infinity breaks the 
remaining spin-boost invariance. 

In the Teukolsky formalism, one considers only space- 
times describing perturbations of a known Kerr metric. 
The physical metric, consisting of both background and 
perturbation, generally does not have repeated principal 
null directions, and does not admit a preferred Kinnser- 
ley tetrad. However, it is nonetheless natural to restrict 
attention to null tetrads of this metric which differ from 
the background's Kinnersley tetrad only at first order in 
perturbation theory. This is Teukolsky's approach. In 
such a perturbed Kinnersley tetrad, all of the curvature 
components except \E'2 vanish at leading order, though all 
may acquire first-order corrections. As a result, all com- 
ponents except ^'2 are invariant under first-order diffeo- 
morphisms, the gauge transformations of the linearized 
theory. Moreover, the residual, first-order ambiguity in 
the tetrad is generated by null rotations with parame- 
ters a and b which vanish at leading order, and spin- 
boosts with parameters c equal to unity at leading or- 
der. The exchange operation, being finite, is disallowed. 
These residual transformations certainly modify the per- 
turbed tetrad, but do not affect the values of the cur- 
vature components ^'o, ^"2 and '^4. Thus, \J'o and ^4 
evaluated in any such perturbed tetrad are invariant in 
both senses, and acquire physical meanings. Specifically, 
using an analysis of the geodesic deviation equation due 
to Szekeres [ll|, [lj| , they are associated with transverse 
gravitational radiation propagating along the null direc- 
tions n° and respectively. 

The components 4'i and in a perturbed Kinnersley 
tetrad are invariant under infinitesimal diffeomorphisms, 
but may be given any value whatsoever by performing 
suitable first-order null rotations on the tetrad. They 
are therefore pure gauge in the Teukolsky formalism. In 
a variety of situations, including both the reconstruction 
of metric perturbations in the first-order theory [T^ and 
a second-order perturbation analysis |16| | , it is convenient 
to use the remaining null rotation freedom to set 

«'i=0 = *3. (7) 

This can always be done, and exhausts the first-order null 
rotation freedom in the tetrad. There is no similar con- 
venient criterion to eliminate the first-order spin-boosts, 
but none is needed in the Teukolsky formalism. 

Let us now focus on the subtleties inherent in apply- 
ing Teukolsky's techniques in numerical relativity. Black 
holes in general relativity are known to be stable . 
Thus, when a numerical evolution produces a space-time 
containing a single black hole, one should expect it to set- 
tle down over time, approaching a quiescent final state. 
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A description via perturbation theory should accordingly 
become increasingly viable over time, although in prac- 
tice it may be far from clear exactly how a given physical 
metric should be so described. One natural response to 
this difficulty would be to search among all possible Kerr 
metrics on the given space-time for one in which the per- 
turbation, the difference between the physical metric and 
that particular Kerr metric, is smallest. However, even if 
it were clear how to conduct such a search systematically, 
it would likely be an extremely difficult procedure. We 
take a different approach here. The key element of the 
Teukolsky formalism provided by the background Kerr 
metric is its Kinnersley tetrad. Rather than search for 
an approximate Kerr metric on space-time whose Kin- 
nersley tetrad will allow us to exploit the Teukolsky for- 
malism, we search directly for an approximation to that 
Kinnersley tetrad. 

Given only the physical metric on a perturbed Kerr 
space-time, a natural way to begin searching for its per- 
turbed Kinnersley tetrad is to seek tetrads for the phys- 
ical metric satisfying Eq. ([TJl. That is, to find a tetrad 
with and vanishing at leading order, simply de- 
mand they vanish to all orders. This eliminates only 
gauge degrees of freedom in first-order perturbation the- 
ory, and there is no obvious problem of principle seeking 
such tetrads at an arbitrary point of a general space-time. 
The components and ^'4 in such a frame cannot gen- 
erally be set to zero simultaneously by further refining 
the choice of tetrad. Even in perturbation theory, these 
components are physical and provide a rough, quanti- 
tative measure of how "far" such a space-time is from 
Kerr |^ ^| . Motivated by their association with trans- 
verse radiation in perturbation theory, we refer to tetrads 
satisfying Eq. (|7|) as transverse null tetrads. We empha- 
size, however, that in regions of space-time supporting 
strongly self-gravitating fields, there may be no precise 
relation between these tetrads and any proper notion of 
gravitational radiation. 

Even in the comparatively simple case of a perturbed 
Kerr space-time, there may be many different transverse 
tetrads. These may differ from one another in two ways. 
First, the transversality condition of Eq. is invari- 
ant under the spin-boost and exchange transformations 
of Eq. (@J whence, whenever one transverse tetrad ex- 
ists, so do infinitely many others. In addition to this 
generic multiplicity, however, there could exist distinct 
transverse tetrads on a given space-time related by non- 
trivial null rotations. This second ambiguity raises a 
more serious problem in our context since even the di- 
rections of the real null vectors in a transverse tetrad 
on a perturbed Kerr space-time may have nothing to do 
with the underlying Kinnersley tetrad. The current pa- 
per addresses this second difficulty. It constructs, up to 
spin-boost and exchange, a unique tetrad on a broad class 
of space-times comprising a non-perturbative neighbor- 
hood of the Kerr geometry in the space of solutions to 
Einstein's equations. We refer to the product of this con- 
struction as the quasi-Kinnersley frame. Mathematically, 



it is an equivalence class {£"•, n", m°, m"} of null tetrads 
which can be transformed into one another using only 
spin-boost and exchange transformations. The set of all 
tetrads at a point is a disjoint union of such equivalence 
classes, and we refer to these generally as null frames. 
Although this nomenclature may be slightly unfortunate 
— the terms "null tetrad" and "null frame" are often 
used interchangeably in the literature — we will adhere 
strictly to it here since the notion of a frame is essential 
to our argument. 

III. THE QUASI-KINNERSLEY FRAME 

A null frame is completely determined by the com- 
mon directions, {i"^} and {n"}, of the real null vectors 
in its constituent tetrads. These real null vectors span 
a time-like 2-plane in the tangent space at each point of 
space-time, while the complex vectors span the orthogo- 
nal space-like 2-plane. Spin-boost and exchange opera- 
tions may change the bases defined by the tetrad on these 
two sub-spaces, but not the sub-spaces themselves. Since 
knowledge of either orthogonal sub-space determines the 
other, the real null directions suffice to determine the 
frame. Mathematically, we can exploit this interrelation 
by associating a 2-form'^ 

Safe = -^[a "-;,]- "If^mf,]. (8) 

to any given null tetrad ,n°- ,m°- ,fh"'). This complex 
2-form is invariant under spin-boosts and only changes 
sign under exchange. In addition, it is self-dual* and 
unit: 

:= i e,b,d E^-^ = +1 S.fa and E^bS^'^-l. (9) 

In fact, all unit self-dual 2-forms arise in this way. For 
each, one can find an oriented null tetrad of vectors 
{1°-, n°- ,m°- ,fh°-) such that S^b takes the form of Eq. 
For a given Y,ab, we have 

V ^ -b \- (10) 

That is, £° and n°- are the unique real null eigenvectors 
of Sab I while to" and to° are its unique complex null 
eigenvectors orthogonal to both and Every unit, 
self-dual 2-form has such eigenvectors, each of which is 



^ Here, and throughout this paper, our convention is that round 
(square) brackets on indices denote the {anti-)symmetric part of 
a tensor. That is, Tat = Tf^ab) + T^ab] ■ 

Although self-dual 2-forms are common in the mathematical rel- 
ativity literature, they are seldom used in numerical work. A 
summary o f th e relevant mathematics may be found in Chapters 
3 and 4 of Il3l . Note, however, that the notation of this paper 
differs in a few details. We have therefore endeavored to make 
the presentation here as self-contained as possible. 
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determined only up to scaling, whence the tetrad is not 
unique. Rather, there is a canonical one-to-one corre- 
spondence between null frames and non-null directions 
in the complex, three-dimensional space of self-dual 2- 
forms at each point of space-time. 

A non-null, self-dual 2-form Eaf, is associated to a 
transverse null frame if and only if 

Cah"'' Scd = -4*J Sab, (11) 

where Cab'^'^ denotes the Weyl curvature tensor of space- 
time. For most T^ab, the right side would include addi- 
tional terms of the form 4^i fh[a rih] and 4\['3 £[q mf,]. It 
is straightforward to show that whenever Scd is self-dual 
on the left side of Eq. , the right side is self-dual as 
well. Thus, the search for transverse frames reduces to an 
eigenvalue equation for the Weyl tensor, viewed as a map 
from the space of self-dual 2- forms to itself. Notably, the 
eigenvalue for a given eigen-form of the Weyl tensor is di- 
rectly connected to the Newman-Penrose scalars in the 
associated null frame: ^'J denotes the ^'2 scalar, which is 
invariant under both spin-boost and exchange, evaluated 
in a transverse frame. 

The above result establishes a one-to-one correspon- 
dence between transverse frames {^°, n", to°, m"} and 
non-null, self-dual eigen-forms T,ab of the Weyl tensor. 
The number and nature of the these eigen-forms are very 
well understood and depend solely on the algebraic 
(Petrov) class of space-time. In particular, there are ex- 
actly three transverse frames at each point of a alge- 
braically general (type I) space-time. Meanwhile, the 
Weyl tensors of type D space-times, such as the Kerr 
solution, have one degenerate eigen-space, leading to in- 
finitely many transverse frames, as well as a distinguished 
non-degenerate eigen-form associated with the Kinners- 
ley frame. Type II geometries admit only one transverse 
frame. In each of these three cases, therefore, the possi- 
ble ambiguity in the transverse frame of interest is finite. 
One can recognize the Kinnersley frame of a type D or 
type II geometry immediately, and the only ambiguity 
occurs in the general case, where there are only three 
possible choices. The most specialized algebraic classes, 
types N and III (as well as type 0, where the Weyl tensor 
vanishes), are very rare and correspond to highly spe- 
cialized fields which are unlikely to arise in numerical 
relativity unless sought specifically. We do not consider 
them further here. 

The central open issue is this. At a typical point in 
the exterior of a black hole formed in a numerical sim- 
ulation, the Weyl tensor will most likely be of type I. 
At late times, perturbation theory is expected to become 
sufficiently accurate to be meaningfully applied and the 
perturbed Kinnersley tetrad, when fixed to first order 
according to Eq. Q, will then belong to a transverse 
frame of the physical metric. In order to implement the 
Teukolsky formalism in this situation, one must discover 
to which of the three transverse frames of the physical ge- 
ometry this tetrad belongs. There are a number of ways 
one might do this in practice but all are based on 




FIG. 1: The Riemann surface Zs is a triple-cover of the Rie- 
mann sphere. The values of Z{S) on its three sheets give the 
three values of tp. The labels in the figure indicate contiguous 
regions of Zs pictured on different sheets. 



the same basic idea. If it is possible to distinguish the 
eigen- value of the Weyl tensor associated with the quasi- 
Kinnersley frame, then its eigen-/orm can be singled out 
using Eq. (|11|) . The frame may then be reconstructed 
from the form using Eq. (|10(l . 

The three eigenvalues of the Weyl tensor at a given 
space-time point are independent of any basis or coordi- 
nate system one might use to evaluate them; they are al- 
gebraic scalar curvature invariants of the Weyl tensor. As 
such, they should be calculable directly from space-time 
geometry. To realize this point explicitly, we calculate 
the characteristic polynomial 

P := - i/V^+ iJ = 0. (12) 

of the linear map defined by the Weyl tensor from the 
space of self-dual 2-fornis to itself. This calculation is 
described in some detail in an Appendix. The roots "0 of 
this polynomial are the three values of in the three 
transverse frames at a generic point of space-time, which 
are related to the eigenvalues by Eq. Hll() . The quan- 
tities / and J appearing here are the well-known scalar 
curvature invariants defined by 

r 1 ( cd ab ■ cd i^r^ ab\ 

-< •— Yg \yab ^cd ~l'~^ab '^cd ) 

T 1 //^ cd cf ab ■ cd ef ab\ ^ ' 

•— gg ((-^ab C/cd ■' L-e/ — iL^ab ^cd ' '^ef ) , 

where *Cab'^'^ '■— \ ^ab""^ Cmn'^'^ denotes the space-time 
dual of the Weyl tensor. These invariants are manifestly 
calculable directly from space-time geometry, with no 
mention of Newman-Penrose tetrads, and show implic- 
itly that the three ^'J are indeed curvature invariants. 
This result may be somewhat surprising at first since the 
curvature component ^2 generally depends entirely on 
the frame used to evaluate it. However, we seek here 
its values only in specific, geometrically privileged, trans- 
verse frames. Since these frames are defined invariantly, 
so are the values of in them. 
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The cubic polynomial of Eq. (|12|l is simple enough to 
solve analytically, but its solutions will involve multiple- 
valued functions of the complex variables / and J . To 
minimize the subtleties raised by these functions, it is 
convenient to express the solution in terms of the Baker- 
Campanelh speciality index |17] S :— 277^//"^, yielding 



(14) 



where W{S) := VS — \/5~~T. The multiple values of 
tjj then arise from the various branches of Z(S). This 
function has branch points of order two at 5 = and S = 
1 arising from the square roots in W{S) and, since W{S) 
does not vanish for any finite value of S, only one branch 
point of order three at 5 = oo. A careful analysis shows 
the Riemann surface Zs for Z{S) has the three-sheeted 
structure depicted in Fig. ^ The three eigenvalues of the 
Weyl tensor arise from the values of Z{S) in these three 
sheets. 

The speciality index was introduced to provide a quan- 
titative measure of whether a given physical metric 
should admit an accurate description in perturbation the- 
ory [T^ . It equals unity if and only if space-time is either 
of type D or of type II. (It ill-defined for space-times of 
types N and III, where I = = J, but we do not con- 
sider these exotic cases here. They could be incorporated 
consistently, however, by working with a multivariate ex- 
pression for Z in terms of / and J.) One therefore expects 
iS ^ 1 at late times in a numerical space-time describing 
a black hole approaching equilibrium. Only one branch of 
Z{S), associated with the leftmost leaf of Zs in Fig. ^ 
is analytic in this limit. Let ip'^ denote the value of ip 
computed in this branch and expand about 5 = 1 to find 



J 



-[-3-f|(5-l)-f...] 



(15) 



Because Z{S) is analytic at 5 = 1 in this leaf of Zg, 
we find a proper Taylor expansion. On the other hand, 
5 = 1 is a branch point of order two in the other two 
leaves, and the expansion there should accordingly be 
done in powers of — 1 rather than 5 — 1. Denoting 
the values of tp computed in these branches by ■0*, we 
expand to find 
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(5-1)1/2 
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(16) 



These two expansions are completely independent of how 
the branch lines for Zg are drawn in the complex 5-plane. 
The branch point at 5 = 1 makes it impossible to distin- 
guish the roots ijj^ from one another in any neighborhood 
of that point, and changing the ambiguous sign of \/S — 1 
indeed merely interchanges these two roots. 

Since the eigenvalues are equal when 5=1, one ex- 
pects they must be associated with the degenerate eigen- 
space of the Weyl tensor on a Kerr space-time, meaning 
tp'^ must be associated with the Kinnersley frame. One 



can check this explicitly by writing the / and J in terms 
of Newman-Penrose components: 



J = ^'o^'2^'4 - *2 + 2^i^'2^3 - ^"0*3 " ^l^i- 



(17) 



In the Kinnersley frame on a type D space-time, only 
^'2 is non-vanishing, and these formulae show that = 
—3 J/ J in that frame. This, of course, is exactly the value 
of ■0° when 5=1. However, the expansion of Eq. p5|l 
remains valid within some finite radius of convergence in 
the complex 5-plane. Within this radius, we can rigor- 
ously identify a unique eigen-form of the Weyl tensor 
satisfying 



Ccd x^O 
ab ^ 



cd 



(18) 



We then define the quasi-Kinnersley frame as that associ- 
ated to S"^ by Eq. (|10|l . This definition is valid whenever 
5 is sufficiently close to unity. 

This quasi-Kinnersley frame may be selected in a num- 
ber of alternative, but equivalent, ways. Definition 3 
of the companion paper [Tlj offers one such alterna- 
tive, essentially using the distinct limits -0° —3 J/ 1 
and t/)* 3J/2I as 5 1 to distinguish the desired 
frame. That procedure will of course select the same 
quasi-Kinnersley frame described above. However, we 
have formulated the definition given here to emphasize a 
point. In numerical work, one has at one's disposal only 
the given physical metric and cannot actually take the 
limit 5 — *■ 1 at a particular point of the numerical space- 
time. Nonetheless, for a given 5 7^ 1, Eq. IjlSI) gives an 
unequivocal value for ■0'^ which may be used to select 
the quasi-Kinnersley frame via Eq. I|18|l . This branch, 
and therefore the associated frame, does not exist only 
in the limit 5^1, but rather in a finite neighborhood 
of 5 = 1 in Zg given by the radius of convergence of the 
series of Eq. ((T^ . Thus, although the quasi-Kinnersley 
frame is selected fundamentally by its limiting behavior, 
it is defined at points of space-time where 5 differs quite 
substantially from unity. To understand the limitations 
of the definition, therefore, we need to analyze the series 
of Eq. (|15|l . Fortunately, the essential features of the se- 
ries are evident in the structure of the Riemann surface 
Zs- They are brought out by the following, relatively 
practical, considerations. 

In practice, the computer cannot calculate ip'^ by sum- 
ming the (infinite) series whose leading terms are given 
in Eq. (|15|l . This series can be very slow to converge 
at points. Rather, it is better to calculate -0° non- 
perturbatively by picking particular branch cuts and 
branches for the square- and cube-roots in Eq. H14() such 
that the series expansion of tp with those choices yields 
■0*^ near 5 = 1. It turns out the principal branches of 
all three radicals, where \/T = 1 = with the conven- 
tional branch lines along the negative real axis yield -0" 
everywhere in the disk |5— 1| < 1. The ambiguous sign of 
\/S — 1 causes no problem in this region because chang- 
ing it, given the branches of the other radicals, merely 
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interchanges the two cube roots in Z{S). Thus, ^/^^ is 
continuous across that part of the real axis joining 5 = 
and 5 = 1. This is a crude reflection of the absence of 
any branch line going to 5 = 1 in the corresponding leaf 
ofZs. 

When |5 — 1| < 1, we have given a simple algorithm to 
calculate xj:^, and therefore the quasi-Kinnersley frame, 
non-perturbatively. Every perturbed Kerr space-time 
studied by the Teukolsky formalism will certainly satisfy 
this inequality everywhere outside the black hole. On 
such space-times, the quasi-Kinnersley frame provides a 
gauge-fixed version of the perturbed Kinnersley frame 
which is manifestly derived solely from the physical met- 
ric. However, other space-times, which in no way repre- 
sent small perturbations of the Kerr geometry, also ad- 
mit a quasi-Kinnersley frame. This structure is therefore 
uniquely and rigorously defined on every space-time in a 
finite neighborhood of the Kerr solution in the space of 
solutions to the Einstein equations. 

It may be possible to extend the definition of the quasi- 
Kinnersley frame to even more general regions of space- 
time using a continuity argument. The critical ingredient 
of the definition which fails when |5 — 1| > 1 is that it 
becomes impossible to distinguish tp^ from tp^ because 
of the branch point at 5 = in Zg. Concretely, if one 
continued to use the principal branches of the radicals in 
Eq. 114|1 to define i/'^ in this broader region, this eigen- 
value, and therefore the distribution of associated frames, 
will be discontinuous at any point of a given curve in 
space-time where S is real and negative. This discon- 
tinuity, however, stems entirely from the branch line in 
Zg. If we instead define a function ^ along such a curve 
which equals "0" up to the point on the curve where S is 
real and negative, and thereafter equals one of the other 
roots, either ip~^ or ip~ , then -0 will allow a continuous 
extension of the quasi-Kinnersley frame along the curve. 
The roots can be found using the principal branches 
of the two square roots in Eq. (|14() . and the branch of 
the cube root where s/T = (— 1 ± i^/3)/2. In principle, 
it appears this extension of the quasi-Kinnersley frame 
to strong-field regions can be carried out coherently (i.e., 
such that the continuous propagation of the frame from 
one space-time point to another is independent of the 
curve used to join them) provided there is no 3- volume 
in space-time where the curvature is everywhere either 
type N or III. In practice, however, if one is interested 
in harnessing the conceptual power of the Teukolsky for- 
malism in numerical relativity, this may not be worth- 
while. This "root-substitution" scheme could incur a 
high computational expense since it can only be done 
retroactively, once the late-time evolution is known, and 
perturbation theory would likely be unreliable anyway 
when |iS — 1| > 1. Nonetheless, there may be situations 
where the numerical overhead could be justified. For ex- 
ample, has advanced a definition of a radiation scalar 
based on exactly this sort of approach to defining the 
quasi-Kinnersley frame in strong-field regions of space- 
time. That argument yields a scalar invariant ^, defined 



throughout space-time, which has a clear interpretation 
in Teukolsky- type radiation zones. 

IV. INITIAL- VALUE FORMULATION 

The previous section has defined the quasi-Kinnersley 
frame from a space-time perspective. However, this lan- 
guage may not be the most convenient for numerical rela- 
tivity. Few numerical codes currently implement a space- 
time covariant form of the Einstein equations, preferring 
a space-|-time split where gravity is encoded in a set of 
fields evolving relative to a fiducial foliation by space- like 
hypersurfaces. These fields include at least the intrinsic 
spatial metric hab and extrinsic curvature Kab of each 
leaf of the foliation. This section describes a simple way 
to calculate the quasi-Kinnersley frame using only these 
initial-value data, thereby offering a direct path toward 
numerical implementation. 

Let f° denote the unit future-directed normal to a 
given space-like hypersurface in space-time. Following 
the convention used in Maxwell theory, a general unit 
self-dual 2-form "Eab can be decomposed into electric and 
magnetic parts as^ 

"^ab ^[aU] + ^'^abcO-". (19) 

Here, eabc '■= tmabc denotes the intrinsic spatial vol- 
ume element and a" is a complex spatial vector with real, 
unit norm. Since Sab is self-dual, its electric part (tech- 
nically, —a"-' /2) also determines its magnetic part, and 
therefore Safe itself, completely. Thus, at each point on 
a given spatial hypersurface, there is a one-to-one corre- 
spondence between complex spatial unit vectors (up to 
sign) and null frames on space-time at that point. 

An equation which determines the (t° corresponding 
to transverse frames follows by taking the electric part 
of both sides of Eq. l(TT|l : 

C%a^ := (£;'^b-iB%)<T'' = 2*Ja". (20) 

This is once again an eigenvalue problem whose solu- 
tions are proportional to the \l/2 scalar evaluated in a 
transverse frame. The spatial tensor Cab featured here is 
symmetric and trace-free, though generally complex. Its 
real and imaginary parts are the electric and magnetic 
components of the Weyl tensor itself, defined relative to 
the given spatial slice by 

^ab • ' ' ^ambn 

21 

^ab • — ' ' ^ambn • — 2 *^amij ^ bn^ 

respectively. Both of these tensors are obviously real, 
and the symmetries of the Weyl tensor imply each is sep- 
arately symmetric and trace-free. 



Fol lowing the convention of Penrose's abstract index notation 
, latin indices are used to denote both space-time and spatial 
tensors. The intended meaning should be clear in context. 
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The Gauss-Codazzi-Mainardi relations of differential 
geometry yield expressions for Eab and Bab in terms of 
initial- value data. The result, which applies regardless 
of whether the space-time metric satisfies the Einstein 
equations, suggests we consider the quantities 



lb 



^Rab 



KK 



b - 
m -^^nb 



K„ m. Kh 



1 uni l-,n 4f^ 

2 "a 



(22) 



Here, K :— Kab is the traced extrinsic curvature, 
while ^Rab is the Ricci curvature of the torsion-free con- 
nection Da compatible with the spatial metric hab- The 
last term in the expression for Sab is the spatial projec- 
tion of the space-time Einstein tensor. When the Ein- 
stein evolution equations hold, therefore, this term can 
be replaced with the stress tensor of matter in the spa- 
tial slice, and thereby expressed in terms of Cauchy data 
for the matter fields. Thus, £ab and Bab are functions 
only of initial-value data when the equations of motion 
hold. Moreover, £ab is naturally symmetric and its trace- 
free part is Eab, while Bab is naturally trace- free and its 
symmetric part is Bab- Thus, one can compute the ten- 
sor Cab in Eq. H2U|) from initial-value data for both the 
gravitational and matter fields. 

The eigenvalue problem of Eq. is easy enough to 
solve in practice — Cab is, after all, only a 3 x 3 matrix 
once coordinates are chosen — but we must once again 
pick from among its (generally) three solutions the one 
corresponding to the quasi-Kinnersley frame. To do so, 
we first find 



/— |C°fcC^a and J ~ —^C"'bC''cC'^a- 



(23) 



These formulas give the fundamental algebraic invari- 
ants of the space-time Weyl tensor as functions of spatial 
initial-value data. One can then easily compute the spe- 
ciality index and, provided — 1| < 1, calculate V''^ using 
Eq. (|14|) . As before, ip'^ can be found using the principle 
branches of the square- and cube-roots in Eq. (|14|l with 
the branch cuts for all three radical functions along the 
negative real axis. This preferred eigenvalue determines 
a preferred eigenvector (Tq in Eq. H2U|I associated to the 
quasi-Kinnersley frame. 

Once the eigenvector CTq has been found, one might 
like to find the elements of some tetrad in the quasi- 
Kinnersley frame. This, too, can be done using only spa- 
tial, rather than space-time, data. The real null vectors of 
an arbitrary tetrad project into a space-like hypersurface 
to define a pair of real spatial vectors. While the normal- 
izations of these spatial projections naturally vary if one 
performs a spin-boost on the tetrad, they define a pair of 
invariant rays in the tangent space at a point of space- 
time. These rays correspond to a pair of real unit vectors 
A'' and i)" which are generally completely independent of 
one another. In particular, they point in exactly opposite 
directions only when the normal t"" to the spatial slice 
lies in the space-time tangent 2-plane spanned by and 
n". The fiducial time foliation used in numerical relativ- 
ity is overwhelmingly likely be "boosted" relative to the 



2-plane defined by a given null tetrad, and only parallel 
unit vectors are actually disallowed. 

Since A° and determine the directions of the space- 
time vectors and n°-, they also suffice to determine the 
frame associated to ct", and thus (t" itself. The explicit 
relation is 



a'' = (1 - A • !>)-' {X" ie'''"' Xb Vc) 



(24) 



Note that when A" and j)" are interchanged, a"' simply 
changes sign, as it should since Yiab changes sign under 
an exchange operation in the corresponding frame. 

To invert Eq. (|24|l and solve for A" and given ct", 
we separate ct" into real and imaginary parts: 



(25) 



In general, neither x'^ nor ?/" is unit, but the normaliza- 
tion condition for a"' demands they be orthogonal with 
norms satisfying — ||t/||^ = 1. Taking the electric 
part of Eq. ((Tn|l then yields 

^a^^abc^^y^ _^a^^abc^^y^ 

X — — tt; and v — — 



(26) 

Finally, the elements of the most general null tetrad 
(^", n", m°, m") in the frame associated to a" take the 
form 



\c\- 



{t" -f A"^) 



(27) 



IIL ^ , — 

where the complex unit projection of to" is 

/t" = (1 A • i>)-i (A" + !)'' + ie"^" Xb vc). (28) 

These vectors are the result of a spin-boost, with param- 
eter c = |c|e**, acting on a preferred tetrad determined 
by setting 1°" Ta — n"" Ta < and m°- Ta — fh'^ Ta < 0. 

Note that in the degenerate case where A° and are 
indeed anti-parallel, the expression for to° in Eq. H27|) is 
ill-defined. However, the limit as t)" —A" does exist, 
though it does depend on how the limit is taken. Should 
this situation arise in practice, it is easy enough to accom- 
modate. The first two of Eqs. H27|l are non-singular, and 
remain unchanged. One can then choose any real unit 
vector f"' in the spatial 2-plane orthogonal to A° = — i)", 
and take 



^abc 



(29) 



This replaces the third of Eqs. (|27|l to give the remain- 
ing two elements of the desired (oriented) null tetrad on 
space-time. Note that there are two arbitrary parame- 
ters, the unit vector and the spin parameter 9, needed 
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to specify in this degenerate case. They are redun- 
dant. Any rotation of f° in the plane orthogonal to A° 
may be compensated by an adjustment to the spin pa- 
rameter 9] as in the non-degenerate case, the spin-boost 
freedom in the tetrad will remain. The only cost of degen- 
eracy is the loss of the natural "gauge-fixing" of the spin 
freedom implicitly used in Eq. (|27|l . That is, we cannot 
specify = by insisting that f" nia be real and negative 
since, in the degenerate case, it vanishes for all values of 
9. In any event, this amounts only to a slight mathemat- 
ical inconvenience. The residual freedom in the tetrad 
is the same regardless of whether A" — —v"" and, more- 
over, this circumstance is unlikely to be realized exactly 
in numerical relativity. 

This section has now shown explicitly that transverse 
frames can be calculated directly from data on an arbi- 
trary spatial slice of a given space-time. In particular, 
one can find these space-time quantities using the fidu- 
cial foliation underlying any numerical evolution scheme 
based on the familiar Arnowit-Deser-Misner equations 
or their various descendants. 



V. CONCLUSIONS 

The central point this paper has argued is that the 
quasi-Kinnersley frame is naturally defined on a broad 
class of space-times. As long as space-time contains no 
finite volumes of very high algebraic speciality, a "root- 
substitution" scheme can propagate the quasi-Kinnersley 
frame continuously from regions supporting fairly weak 
fields, where there is absolutely no ambiguity in its defini- 
tion, to arbitrary regions supporting much stronger fields. 
To avoid potential subtleties inherent to this scheme, one 
can restrict attention to situations where |iS — 1| < 1 
throughout space. In this restricted case, the quasi- 
Kinnersley frame is particularly easy to find. In the still 
more restrictive limit 5 — *■ 1, the quasi-Kinnersley frame 
always exists, depends only on the physical metric, and 
identifies a "gauge-fixed" version of the perturbed Kin- 
nersley frame satisfying Eq. ((T)). This provides essen- 
tial guidance needed to deploy the Teukolsky formalism 
where it is expected to be relevant. This feature could 
be useful in extracting geometrically-invariant, physically 
meaningful information from the raw numerical data of 
late-time black hole simulations. 

The mathematical structure discussed above could also 
be useful in constructing observer-independent scalar 
functions in more general regions where perturbation the- 
ory may not apply. Although the physical meaning of 
these scalars may be somewhat obscure in regions sup- 
porting self-interacting fields, they can still be used, in 
principle, as gauge-invariant observables and thus help 
clarify invariantly the structure of simulated space-times. 

It does not seem possible to establish concrete, analytic 
limits on the regime where perturbative techniques are 
justified. Their viability must be tested by some other 
means. The internal consistency checks of the Lazarus 



project [l9j, where late-time evolution is done by extract- 
ing Teukolsky quantities from numerical data on a given 
time-slice in a non-linear evolution and thereafter using 
the Teukolsky formalism's much more stable evolution 
equations, provide a natural way to do this. These con- 
sistency checks ensure the physical content of the very 
late-time data does not depend sensitively on when the 
extraction took place. This paper has argued that one 
can begin attempting extractions whenever one first finds 
|5 — 1| < 1 everywhere outside a single black hole. In 
practice, it appears this can happen surprisingly early in 
a full numerical evolution pjj . 

In addition to transverse frames, it is possible to con- 
sider several other geometrically-preferred frames on a 
typical space-time. For example, frames whose real 
null directions are two of the (generally four) principal 
null directions of the physical metric are quite natu- 
ral to consider. The connection between such principle 
null frames and the transverse frames considered here is 
discussed much more thoroughly in the companion pa- 
per ^3 (see especially section II. F). One might won- 
der whether techniques similar to those outlined here for 
the quasi-Kinnersley frame could be developed to pick 
some of these other preferred frames on a generic ge- 
ometry. The answer, however, is probably in the neg- 
ative. The essential feature which has distinguished 
the quasi-Kinnersley frame in this paper is the analyt- 
icity of ■i/'^ near 5=1. Even the other two transverse 
frames cannot be distinguished from one another near 
algebraic speciality because of the branch cut joining 
their associated leaves of the Riemann surface at this 
point. One can show using standard perturbation theory 
of linear operators that, since "0* fail to vary analyti- 
cally near speciality, their associated frames also fail to 
do so This non-analyticity, and our inability to dis- 
tinguish these frames invariantly, arises precisely because 
we have perturbed away from a situation where the nor- 
mally distinct frames degenerate. All of the other obvious 
geometrically-preferred frames, such as those based on 
principal null directions, also degenerate in algebraically 
special space-times, and thus will share this problem. 
The quasi-Kinnersley frame is unique in that it varies 
analytically when perturbed, and is therefore arguably a 
singularly useful structure in the sort of analyses contem- 
plated here. 

As discussed in the Introduction, a second step must be 
taken to complete the construction begun in this paper 
and allow an application of the full range of Teukolsky's 
techniques to the results of numerical simulations. This 
step will break the spin-boost-invariance of the current 
results and choose a preferred quasi-Kinnersley tetrad 
within the quasi-Kinnersley frame. Unlike the frame- 
fixing procedure described above, which has been entirely 
local, done separately at each space-time point, fixing the 
scalings seems likely to require a global approach. Setting 
e = 0, as for the Kinnersley tetrad on Kerr space-time, 
can still be done generically, and will continue to help 
pare down the possibilities for a quasi-Kinnersley tetrad. 
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However, as in Kinnersley's case, this will not determine 
the tetrad uniquely. In fact, since the e = condition 
ties together the scalings of tetrad elements along the in- 
tegral curves of it would seem that fixing the scaling 
close to the black hole would involve making choices in 
the evolved fields at the outer boundary at much later 
times. Since long-time evolutions of the full, non-linear 
equations of general relativity are currently unavailable, 
it seems likely an indirect approach to this problem will 
be required for the near-term. It may be possible to found 
such a approach on a scheme to approximate additional 
structures of the Kerr geometry underlying a given per- 
turbed metric using only the physical metric itself. Work 
is currently underway to explore whether such an ap- 
proach could yield a practicable solution to the problem 
of finding the quasi-Kinnersley tetrad. 
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APPENDIX A: THE CHARACTERISTIC 
POLYNOMIAL 

The characteristic polynomial of Eq. (|12|l plays a key 
role in identifying the quasi-Kinnersley frame. For def- 
initeness, this Appendix will derive it. The derivation 
follows closely the discussion given in Chapter 4 of [3. 

The space of self-dual 2-forms at a point of space-time 
is a three-dimensional complex vector space. To see this, 
we recall that in Lorentzian signature, any 2-form Qab 
may be written as a sum of its self-dual and anti-self- 
dual parts: 

nab = ^^ab + ^^ab, (Al) 

where ^flab '■— \ [^ab Ti*^^a&] are the self-dual (upper 
sign) and anti-self-dual (lower sign) parts of Qab- One 
can check the 2-form 'H^ab {~^ab) is (anti-)self-dual us- 
ing the result **^ab ■— *{*^ab) = — f^afc, which holds 
for the double dual of any 2-form in Lorentzian signa- 
ture. The space of all 2-forms in four dimensions is six 
dimensional, but the decomposition of Eq. (|A1|) splits 



it into two subspaces which are complex conjugates of 
one another. That is, we see by taking the conjugate of 
the definition of self-duality in Eq. that the complex 
conjugate of a self-dual 2-form is anti-self-dual. These 
two subspaces cannot intersect and together form the 
whole of the space of 2-forms, whence each must be three- 
dimensional. 

The contraction of any self-dual 2-form Sab with the 
Weyl tensor will yield another self-dual 2-form. This hap- 
pens because the Weyl tensor has the remarkable prop- 
erty that its duals on either pair of indices are equal to 
one another: 

cd ^* cd 1 cd { \ 0\ 

^ab — ^ ab ■ — 2 <^'' ^mn ■ \^^) 

Thus, the dual of the contraction Cab'^'^ '^cd can be moved 
to act on the second pair of Weyl indices and then to 
act on Scd, producing a factor of and showing the 
contraction to be self-dual. Thus, the Weyl tensor defines 
a linear map from a three-dimensional complex vector 
space to itself. 

In Eq. H19|) . we have shown explicitly how the three- 
dimensional space of self-dual 2-forms Y,ab at a point of 
space-time can be described in terms of complex spatial 
vectors da ■— Sabf'' in an arbitrary spatial hypersur- 
face. These later objects also form a three-dimensional 
complex vector space. Indeed, the map induced by the 
Weyl tensor on self-dual 2-forms is completely equivalent 
to the map C^b on complex spatial vectors defined by 
Eq. H20|) . For conceptual simplicity, then, we may calcu- 
late the characteristic polynomial of the spatial map C°fc 
rather of than the space-time Weyl tensor restricted to 
self-dual 2-forms. 

For the moment, let Cf, be any map on a three- 
dimensional vector space. Using standard techniques in 
linear algebra, one can show the its characteristic poly- 
nomial can be expressed in terms of the traces of C^b 
and its powers in the form 

= P{c) = det [C - c5] 

^c^- (C) c2 + i [(C)2 - c (A3) 
-i[(C)3-3(C)(C2)+2(C3)], 

where S°-b is the usual identity map and (C") denotes 
the trace of the n"^ power of C^b- Returning now to the 
particular case of the Weyl tensor, we see immediately 
that these coefficients simplify considerably because C°fc 
is trace-free. Moreover, Eqs. (|23|l relate the remaining 
traces to the curvature invariants / and J. This reduces 
the characteristic equation for the eigenvalues of the Weyl 
tensor to the form 

- Ic + 2J = 0. (A4) 

Since, by Eq. 12011 . the eigenvalue c is twice 'I' J (in a 
transverse frame), we recover the polynomial of Eq. 112|) 
for the values of tp. 
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